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ABSTRACT 

An iterative method is presented for solving the problem of a uniformly rotating, self- 
gravitating ring without a central body in Newtonian gravity by expanding about the 
thin ring limit. Using this method, a simple formula relating mass to the integrated 
pressure is derived to the leading order for a general equation of state. For polytropes 
with the index n = 1, analytic coefficients of the iterative approach are determined 
up to the third order. Analogous coefficients are computed numerically for other poly- 
tropes. Our solutions are compared with those generated by highly accurate numerical 
methods to test their accuracy. 

Key words: gravitation methods: analytical hydrodynamics equation of state 
- stars: rotation. 



1 INTRODUCTION 



Motiv at ed in pa r t by the r ings of Satur n, Kowalewskvl 
i| 18851 ). iPoincarel l| 18851 ) and iDvsonl (|l89l 1 18931 ) studied 
amongst other things, the problem of an axially symmet- 
ric, homogeneous fluid ring in equilibrium by expanding it 
about the thin ring limit. In particular, Dyson provided a 
solution to fourth order in the parameter a = a/b, where 
a provides a measure for the radius of the cross-section of 
the ring and b the distance of the cross-section's centre of 
mass from the axis of rotation. An important step toward 
un derstanding ring s with o ther equations of state was taken 
bv lOstrikeil l|l964al lbl ll965l ). who studied polytropic rings to 
first order in a and found a complete solution to this order 
for an isothermal limit. 

Numerical methods were developed to study 
such rings and t heir co nnection to the M aclaurin 

Sugimoto 



spheroids (IWoni 



Eriguchi 



198ll; 



Eriguc hi fc Hachis 
2003d ). With numerical 



1974 : 

1985; lAnsorg. Kleinwachter fc Meinell 
methods, it was also possi- 
ble to treat the problem of non-homogeneous rings 
and even with in the framework of General R e lativity 
ilHachisul 1 19861: lAnsorg Kleinwachter fc Meinell l2003bl : 
Fisc her. Horatschek fc Ansordl2005l ). 

Through the use of computer algebra, we were able to 
extend Dyson's basic idea and determine the solution to 
the problem of the homoge neous ring up to the order a 20 
|Horatschek fc Petrofj 120081 ) . In this paper, we present an 
iterative method for performing a similar expansion about 
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the thin ring limit for arbitrary equations of state and a 
number of general results are derived confirming and gen- 
eralizin g work that had already been published bv lOstrikerl 
(|1964bl ). The application to polytropes is considered and or- 
dinary differential equations (ODEs) are found that allow 
for the determination of the mass density. A closed-form so- 
lution can only be found if the value of the polytropic index 
is n = 1, and such rings are considered to the order a 3 . For 
other polytropic indices, the ODEs are solved numerically so 
that results from the approximate scheme can be compared 
to highly accurate numerical results for a variety of equa- 
tions of state. The numerical solutions considered here are 
taken from a mult i-domain spectral program, muc h like th e 
one described in lAnsorg. Kleinwachter fc Meinell (|2003al ) . 
but tailored to Newtonian bodies with toroidal topologies 
(see lAnsorg fc Petrofj 120051 for more information). The so- 
lutions obtained by these numerical methods are extremely 
accurate and thus provide us with a means of testing the 
accuracy of the approximate method. 



2 APPROXIMATION SCHEME 

Numerical results indicate that, independent of the equa- 
tion of state, the shape of the cross-section of a uniformly 
rotating ring tends to that of a circle in the thin ring limit, 
i.e. the limit in which the ratio of the inner radius Qi to the 
outer one q tends to 1. This suggests that, as in the homo- 
geneous case, a Fourier expansion of the quantities involved 
will yield coefficients related to the thinness of the ring. This 
last statement can be made more precise when the explicit 
approximation scheme is presented. 
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Figure 1. A sketch providing the meaning of the coordinates 



Let us begin by introducing the constant b and polar- 
like coordinates (r, x, f) related to the cylindrical coordi- 
nates (g, z, if) by 



g = b - r cos x, 



z = r sin x, 



(1) 



For a given value of tp, constant values of the coordinate r 
are circles centred about (g = b, z — 0) and x measures the 
angle along any such circles. Fig. [1] provides an illustration 
of the coordinates. The surface of the ring in cross-section 
can be described by a function r = r s (x), which we expand 
along with the mass density fx and squared angular velocity 
fi 2 , 



^2^2 Vik(y) cosikx)^ + o{a q ) 

i=0 fc=0 



= 7rG/i c I ^ ilicr 1 + o(a 



9 + 1' 



(X)=» 1 + cos(fc X K + 



and 



where we have introduced the dimensionless radius 

r 

V ~ - 

a 



and the parameter 



(2) 



(3) 



(4) 



(•>) 



(6) 



which tends to in the thin ring limhfl The quantity /i c is 
chosen to be the mass density at the point r = and does 
not represent the density's maximal value, although it will 
not differ significantly from it in general. The absence of sine 
terms in the Fourier expansions is as a result of the symme- 
try with respect to the equatoria l plane, which is k nown to 
hold for stationary solutions (see lLichtensteinlll933r ). 

In this section, we present a method for finding fi q k, 
Qq+i and f3 q k given that the previous terms in a % are known. 

The idea used in Dyson's approximation scheme for ho- 
mogeneous rings makes use of the Poisson integral to de- 
termine the gravitational potential in terms o f the (still un - 
known) function r s along the axis of rotation (|Dvsonlll89d ). 
This is only possible since the mass density is completely 
determined for homogeneous matter once the shape of the 



1 Note that the expansion for r s contains terms with k = 
in contrast to the analogou s expression for homogeneous rings 
( Horatschek & Petroff 2008), where the summation index k runs 
only from 1 to i. The reason behind this will be discussed shortly. 



ring is given. In general, however, it is necessary first to de- 
termine n to the desired order before being able to perform 
the integral. Here we obtain ordinary, second order differen- 
tial equations for n q k(r) by applying the Laplace operator 
to the integrated Euler equation 



U + 



dp' 



M(P') 2 



1 nV = v , 



4nGn + V 2 



dp' 



- 2Q 2 = 0, 



(7) 

(8) 



where Vb is the constant of integration and p the pressure. 
Applied to a function / = f(r,x), the Laplace operator in 
the coordinates (r, x, f) reads 



9V 
dr 2 



Idf 



r 2 dx 2 

1 / d f sin v d f 
r cos x ) ( cos x 



dr 



dx 



Expanding (J5J) in terms of a and requiring that the 
equations be satisfied for each power in a and each term 
in the Fourier expansion then results in ODEs for fi q k(y) 
once an equation of state has been specified. These func- 
tions must be regular at the origin and chosen such that 
Moo(0) = 1 and /iifc(0) = for all other i and k so as to 
be consistent with the choice ^i(0, x) = For k = 0, this 
condition suffices to determine the function uniquely. For 
k — 1,2, ... ,q, the remaining constants in the solution of 
the ODEs are found by requiring that the pressure vanish 
at the surface. Demanding this for each of the coefficients 
in a Fourier expansion, provides q + 1 equations for the re- 
maining q constants. The additional equation can be used to 
determine /3 9 o- It may come as a surprise that a term with 
k — was included in the expansion here for r s , since no 
such term is needed in the homogene ous case. The scale in- 
variance mentioned toward the end oflHoratschc k fc Petrofil 
(2008) means that a need never be specified in that cas^E 
Here, however, the first zero of noo determines the value of 
a, which can then be 'corrected' order for order in a via the 
coefficients /3io. 

Now that the density has been determined and the num- 
ber of unknowns coincides with that of the homogeneous 
case, we can proceed as we did there, a detailed descr iption 
of which can be found in iHoratschek fc Petrofil (120081 '). The 
freedom one has to choose the origin of the coordinate sys- 
tem, r — 0, is used in requiring that it coincide with the 
centre of mass 



2ir r s ( x ) 



J J r 2 /i(r, x) cos x dr dx = 0. 



(9) 



The potential in the vacuum region is determined from the 

2 Polytropes with n = 1, which will be treated shortly, also con- 
tain an interesting invariance. If £/(x), /^(x) = K^/p(x), £7 and 
Vo are solutions to the Poisson and integrated Euler equations, 
then so are aU(x), a/x(x) = K\J a 2 p(x), and aVo, where a 

is an arbitrary scaling factor. This invariance is reflected in the 
fact that the value for a is independent of /i c for polytropes with 
n = 1. 
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Poisson integral via a rather involved procedure and reads 
[7 out = -2-kG^o 2 \ S2 cP^ff-'Atli + o(a q )j , (10) 
where 

dip 



l_d 

6dfe 



z 2 — 2bgcosip 



(11) 



and the terms A; result from an expansion of the Poisson 
integral. 

Making use of expansion formula? for Uout can be 
evaluated along the surface of the ring r s and expanded in 
terms of a. The Euler equation then tells us that 



U ou t(r a ) - if2 2 (a/(j - r s cosx) 2 = Vo 



(12) 



holds. Considering the a q term and equating the coefficients 
of the cos(fcx), k = 1, 2, . . . , q to zero then provides equa- 
tions for determining the remaining unknowns. 

It should be noted that the expansion coeffi- 
cients depend on powers of lncr in general. Because 
lim CT ^o <j(\na) a = for all a, this dependence does not pose 
a problem for the iteration scheme. 



By multiplying (|15p by 7Vfi c a 2 by and integrating from to 
1, one finds that the mass M to leading order can be related 
to the derivative of /loo at the point y = 1 according to 



M = 4tt V 



c a 2 b J nooydy = — na 3 /j, c /cr 



dh 



dy 



(20) 



y=i 



A particularly interesting relation involving the square of 
the mass can be derived by considering the integral over the 
pressure p, which we first expand 



(21) 



P(r,x) = G/4a 2 ^^Pikiy) cos(fcxK + o(a 9 



To leading order, upon taking (|15p into account, the integral 
over p reads 



2tt r a ( x ) 

P :— 2n J J pr(b — r cosx) dr dx, 
o o 



which to leading order is 



P = 4TT 2 Gfl, 



1 

abj 



Pooy dy 



(22) 



3 GENERAL RESULTS TO FIRST ORDER 

The approximation scheme described above allows us to 
draw certain conclusions even without specifying the equa- 
tion of state, thu s generali z ing res ults that were published 
for polytropes by lOstrikerl (|l964bh . To leading order in a, 
where nothing depends on the angle \, the ring (here a 
torus) is equivalent to an infinitely long cylinder, a prob - 
lem that was s tudied by IChandrasekhar fc Fermil (| 1953f ) ; 
lOstri kcr (1964aj). Upon introducing the pressure function 



„ 2^ 2 4, f dpoo 2 , 
-2tt G^ c a b J -^-V d y 

o 



-2-k 2 G^a'b J 



d/ioo 2 j 
Voo—. — y dy 
dy 



8-K 3 Gfi 2 a 4 b Hooy [ I MooJ/'dy' I dy 



(23) 



h 



dp' 



(13) 



and expanding it as with \x in © 

/ 1 i 

K r >X) = Gp, c a 2 i^^2h ik (y)cos(kx)(J l + o(a q 
\;=0 fc=o 

equation ((Sj reads 

■j-z + - -J-) h o +4n(j, 00 = 
dy* ydy) 

to the lowest order in a, since the integrated Euler equation 
tells us that 

SI 2 /G[i c = o(a 2 ) 



(14) 



(15) 



must hold, i.e. 



n = fii = o. 



(16) 



(17) 



At the surface of the ring r = r s , the pressure vanishes, 
corresponding to h(r — r s ) = 0, and we thus find 



and 



/3n = -hi 



ftoo(l) = 
/ dhoo 



(18) 
(19) 



4-K 3 Gfi 2 a 4 b / nooydy 



GM 2 



Airb 



Numerical examples demonstrating how AivbP/GM 2 ap- 
proaches 1 in the thin ring limit for various equations of 
state can be found in Fig. H 

The terms from the expansion (|10[) of the potential in 
the vacuum that play a role up to first order are 



Ai = 2 Hooydy 



Ma 
2n 2 fj, c a 3 



and 



A 2 = 



gMa 2 

27T 2 /i c a 3 



TY 2 fi c a 3 I 
9 := ~nila~ I dV - 



(24) 



(25) 



The coefficient fi 2 from the expansion of the square of the 
angular velocity is 



with 



y=l 



2 = Ai(l + A - 2/3n) + 2A 2 /cr 



ln(8/ff) 



(26) 



(27) 
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Figure 2. Numerical examples demonstrating how inbP/GM 2 
tends to 1 in the thin ring limit for various equations of state, cf. 
equation (123 b - 



and the constant of integration from the Euler equation is 

1. 



Vo 



GM /5A + 9 
2tt6 



+ 9-Pi 



(28) 



The term g — (3u appearing in the above equation can be 
treated further by considering the rotational energy T and 
potential energy W and making use of the virial identity 



= 3P + 2T + W 

„ GM 2 ( „ 1 
= 3P+ — ^-2 



27T 2 G^i 2 a 4 b I [loohaoydy. 



(29) 



By restricting ourselves to the polytropic equation of state 
(see (|34[) ). we can rewrite the above integral to read 



r 

/ 



Hoohooydy - 



(n + 1)A> 



l/n-l 



Ga 2 



1 

= (n + l) J 



Pooydy = 



/l+l/n j 
Moo v d v 



(n + l)M 2 



(30) 



167T 3 ^ a 4 & 2 > 



where the last step follows from ([23} ■ Putting this expression 
into (|29l) and using (|23[) again then yields 



= 



8 



(31) 



Taking into account 1 — Qi/ Qo = 2a, which holds to leading 
order, we can use (|26[) to write 



,3o2 



+ ln 1 



GM 

and (|28l) can be written as 

AnbVo ,_ ^ y 



5GA/ 



In [ 1 - i» 



5 — n 
20 



+ lnl6 



- In 16 



(32) 



(33) 



for poly tropes in the thin ring limit. Similar equations can 
be derived for J (angula r momentum), T, P and via the 
virial identity for W fsee lOst rikcr 1964lJ). These equations 



5GM 




Figure 3. Numerical ring sequences (solid lines) for homogeneous 
matter (upper curve) and polytropes with n = 5 (lower curve) are 
plotted for rings approaching the thin ring limit. The asymptotic 
behaviour as given by equation l|33| l is indicated by the dashed 
lines. 



also hold for homogeneous b odies (n — 0), as was shown in 
Irloratschek Sz PetrofJ (|2008l ). Numerical examples demon- 
strating the behaviour (|33[) are provided in Fig. [3] 



4 MASS DENSITY FOR POLYTROPES AT 
THE ZEROTH ORDER 

The polytropic equation of state is 

p = Kn 1+1/n . (34) 

For large/small polytropic indices n, the equation is referred 
to as 'soft'/'stiff ' and as n tends to zero, /j, tends to a con- 
stant. From now on, we shall use the terms 'homogeneous 
matter' and 'n = 0' interchangeably. For polytropes, (|8]) be- 
comes 



4irGfJ. + K(n + 1)V 2 (V 



- 2£T = 0. 



(35) 



Instead of our coordinate y, we are now going to make use 
of a new dimensionless radial coordinate, applicable to poly- 
tropes 



G2 Mc 



To lowest order in a, and upon introducing 

l/n 



(j, := (j, 



and the expansion 



fi = At, 



l/n 



Vi = fc = 

equation (|35[) reads, cf. p5|) . 



^1 

da; 2 



+ 



1 d 

x dx 



Moo + 



47T 

n + 1 



~ n 

Moo 



0. 



(36) 
(37) 

(38) 
(39) 



This equation is sometimes referred to as one of the general- 
ized Lane-Emden equations (of the first kind) and solutions 
to it h ave been derived and studied in e.g. lGoenner fc Havaj 
(2000). No solutions other than for n — 1 have been found 
for our particular parameters in closed-form and a discus- 
sion using symmetry tran sformations suggests that they do 
not exist, (Goenner 2001). We thus concentrate in the next 
section on the special case n = 1. 
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5 ANALYTIC SOLUTION FOR POLYTROPES 
WITH n = 1 

5.1 The Zeroth Order: ct° 



We rewrite (I39|l for n = 1, remembering that now fl = [i, 

Moo + 2nn Q0 = (40) 



d 2 + l_d_ 

Ax 2 x dx 



and can immediately write down the general solution 



Moo = Ci Jo(V2ttx) + C 2 y (V27ra : ) 



(41) 



where J n is a Bessel function (of the first kind) and Y n a 
Neumann functio n (also called a Bessel function of the sec - 
ond kind), see e.g. lPrudnikov, Brvchkov fc Marichevl (| 1990T ) . 
The condition fi(r = 0, \) = Mc tells us that Ci = 1 and 
C2 = 0. The first positive zero of Jo determines value for a 
from @. We refer to the fcth positive zero of the nth Bessel 
function as j„k and can then write 

ax/r =: a = joi/V^ = 0.959... (42) 



5.2 The First Order: a 1 

The unknown quantities we have to solve for are fj,io(x), 
Hn(x), /3io, /3n, and S72 - From (|8]), one finds the differential 
equations 



d 2 Id. 



(43) 



and 



d , 1 d "\ , fn M 1 d ^oo 

— + r— )Mii + (27T--)u 11 = -^. (44) 



Ax 2 x Ax 



Considering only solutions that vanish at the point x — 0, 
so as to maintain our choice fi(0) = pi c , we find 



fiio = 



and 



Mil = C3J1 + — — (VZnxJo - Ji) , 
2lm V / 



(45) 



(46) 



2.7 

where the argument of the Bessel function is always \ /r 2irx 
unless otherwise specified. The requirement that the density 
vanish at the surface of the rings determines 

P10 = (47) 

and relates the constant C3 to the surface function 

l + 2jo 2 i/3ii 



C 3 



2joi 



(48) 



The constant fhi is determined by stipulating that the cen- 
tre of mass coincide with the point (q — b, z — 0) as in 



n _ 4-jQ 2 l 

011 -M- 



(49) 



Recalling the definition A := In — — 2 one finally obtains 



«2 



2Ji(joi)(A+l) 
joi 



(50) 



from (HJ). 



5.3 The Second Order: ct 2 

To second order, the unknown quantities that have to be 
solved for are (J,2o(x), ^2i(x), "22(2), P20, f%i, 022, and f2 3 . 
The ODEs describing the mass density now read 

(51) 



+ 



1 / dfii 



2a \ Ax 



IJ.11 

x 



+ 



x dupo 
2a 2 da; ' 



— - — 

Ax 2 x Ax . 



and 



d 2 Id. 



+ 



1 / dm 



x 



+ 



M21 



~~9 M22 



x Ajio 



1 Auiq 
a Ax 



= (52) 



(53) 



2a \ Ax x J 2a 2 Ax 
The solutions vanishing at x = are 



^ 2 n nx 1 
M20 = — (1 - Jo) + - 



3tyx' 



3m 



Jo + 



V2tyx ( -To - 



Jdi 



M21 
and 



— C4J1 



M22 = C5 J2 + 



+ 



4 

TlX 



5 + 37ra; 2 
foi 2j'oi 



■h 



x2 \k 2" 1/1 



(54) 
(55) 



(56) 



The constants C4 and C5 can be related to the surface 
function by requiring that /i( r s) = hold independently for 
the coefficients in front of cosx and cos2x- The result is 



C*4 = JOl/fel 



and 



C5 =2 



•2 Q 

• 2 a J 01 I 



P20 — 3 1 ^ 

J01 J01 



(57) 



■ (58) 



11 

W01 



Requiring the same of the coefficient in front of cos Ox gives 

Evaluating © tells us that 

foi = =^ M21 = 0. (60) 
The values for the remaining constants follow from (|12[) : 

fis = (61) 

and 



02 



1 5(A + 3) 1 
4io 4 i 2j 2 x 64' 



(62) 



5.4 The Third Order: cr d 

The third order is the final one to be presented here, but 
the iterative scheme can be applied up to arbitrary order 
assuming that one is able to solve the differential equations 
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for the mass density and perform the necessary integrals. 
The ODEs that result for this order are 



d 2 Id. 

j — 9 H — ^30+ = 0, 

dar x dx , 



(63) 



— - — 

dx 2 x dx , ^ ' 



27T J fl3 



3 / _d_ \ 1 

4a 2 \ di J a 



dfi2 



dx 



3x 2 d/ipo 
4a 3 da; 

-— + - 

2 da; x 



l>22 



2 1 d \ / 2 . 

+ -— I M32 + 2 7T - — U 32 = 



\ dx 2 a; da; 



(64) 
(65) 



and 



1 d 



da; 2 x dx , 



271 " ~ /^33 = — 



a; 2 dpoo 



4a 3 da; 



1 / d \ 1/ld 1", 

+ te 2 "l a; da:- 1 J^ 11 + Sl2d^-x .^ 22 - 



(66) 



The solutions to these equations obeying the requirement 
A*?fc(0) = are 



M30 = 0, 



AJ3i = C 6 - 



(67) 



(A + 1) Ji(joi) 9 7T (j 2 i - 8) a; 2 



2j'o 2 i 



64 



Joi 



+ 



+ 



Joi + 16jgi (5A + 4) - 16 
256J 3 ! 

0F(A+1) Ji(joi)x _ 15 V2n 3/2 x 3 
2^(joi + 16jg 1 (5A + 4)-16)a; 



(68) 



256J& 



J 2j 



M32 = C7J2 

and 

vra; 2 (l6j 2 j (5A + 2) - 40tt3; 2 + 272 + jt n ) 



(69) 



M33 = 



512j 3 



Jl 



01 

2 



57r 2 a; 4 7T (jtn + 8j 2 i(10A + 7) + 16) x 2 



(70) 



+ c 8 J 3 



64j 3 i ' 512$! 

3(joi + 16j 2 1 (5A + 3) + 208) 
256$! 



The constants 030, Ce, CV and Cg are determined by 
requiring that [i(r = r B ) = hold independently for the 
coefficients in front of the cos Ox, cosx, cos2x and cos3x 
terms. The result is 



030 = 0, 



2j 



„-3 
JOI 

512 



5 . a 39 A 
14A + 9 40A + 37 



32joi 



16j 3 i 



C 7 = 



(71) 



(72) 



and 



Cs 



Jo 2 i - § 



^ f 5 X JL fl UU 77 V- 3 

Y536 -I 32 A + / 333 + 768 ) Joi 



,35, 11 

- 1 77i A — J 



1 ^10A + 181 ^ 119 

301 



48 7 6$! 



(73) 



Equation (|9} yields 

9 . 31 
031 = — A + 



32 128 $1 



1/13, 5 
-A + 



1 /5 



:A + 9 



and (fT2j) gives 

032 = 0, 

n 4 = -2 



(A + l) ; 
Joi 



Ji(joi) 2 + 



joi V 8 2 



1 /II 17 



32 



(A+l)joi 



(74) 



(75) 



(76) 



and 



(7') 



5.5 Physical Parameters 

The shape of the rings to third order that results from equa- 
tion ((4]) is compared to numerical results of the correspond- 
ing radius ratio in Fig. [4] For thin rings, the numerical and 
third order curves are indistinguishable. As the radius ratio 
is decreased, the numerical results show that the outer edge 
becomes pointier, right up to the mass-shedding limit for 
the value Qx/q = 0.25322 .... For such a ring, a fluid par- 
ticle rotating at the outer rim in the equatorial plane has a 
rotational frequency equal to the Kepler frequency, meaning 
that it is kept in balance by the gravitational and centrifugal 
forces alone - the force arising from the pressure gradient 
vanishes. The shape of the ring with the cusp that forms 
for mass-shedding configurations is not well represented by 
a small number of terms in our Fourier series. 

Using the results of the last subsection, we write down 
expressions for various physical parameters and can use 
them to verify that the virial identity 3P + 2T + W = 
is satisfied to each order in a. For convenience, we first in- 
troduce dimensionless quantities, valid for any polytropic 
index n > 0: 



M 
M 

P 
P 



J 
J 



G 2 fi~ c 



Ki 



n 



T 
T 



W 
W 



- 



(78) 



K2 



where J refers to the angular momentum. Up to and includ- 
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0.25322 



Q Qo 

Figure 4. Meridional cross-sections of polytropic rings (n = 1) 
with varying radius ratio are shown to third order (solid lines) in 
comparison to numerical results (dashed lines) for the same radius 
ratio. At the value Q\/ Q a = 0.25322 . . ., the rings reach a mass- 
shedding limit, as is evident from the numerical cross-section. 



ing third order one finds 



M = 



2nj 01 Ji(Joi 



+ 



2% (j i + 28j t i + 32 ji^ A - 16) Ji(joi) 



01 



(79) 



27Tj l (A+l) Ji(joi) 



, = J^otKA + l) f j| 

301 cr 2 



(A + l)ioiJi(joi) 



1 



[3j 4 i(A + l) 



(80) 



P = 



128 (A + l) 

+ Joi (96A 2 + 324A + 232) - 624A - 664j } , 



joi Ji(joi) 



2 3oi 



+ 



^[-640joi (A+l) Ji(joi) 



(81) 



+ 10 (joi - 8j 2 ! + 128A + 136) ]a\, 



T= ^HjoiV 3oi\ - (A+l) j i 



[ - 2 (A + l) 2 J! (joi) joi + ^ (jo 4 i(A + 1) 
+ 2j ( 2 1 (8A + 9)(2A + 3) - 112A- 132^] a \ 



(82) 



n 2 0.4 - 




Figure 5. The square of the dimensionless angular velocity is 
plotted versus the radius ratio for rings with polytropic index 
n=l. 



and 



W - 



2~;n! + +a) 2 <j - J(2A + 5) 



joi (2A + 5) (A+l) Ji (joi) 



^ (j^ (2A + 5) + 4j ( 2 i (16A 2 + 42A + 21) 



(83) 



+ 160A + 144 



In the derivation of the above expressions for P and W, we 
have made use of the identity 



20 ( ^, | 2, 2, | -j 2 i) = 3j 2 i 2F3 ( 3, 3, 7 -;-j 2 01 



P F q (a 1 ,a 2 , . . . ,a p ; 61, b 2 , ■ ■ ■ , b q ; z) 



for the Gauss hypergeometric function 

) 

(ai) fc - M fc ■ ■ ■ M fc z k 

; __ (M*-(M fc ---(M fe k\ 

with the Pochhammer bracket 



E 



(84) 



(85) 



(a) k :=a(a + i)---(a + k-i), 



(0)0 : = !• 

A proof of (|84l) can be found in Appendix [A] 

In order to gauge the accuracy of the expressions listed 
above, some of them are plotted to first and third order in 
comparison to numerical values in Figs [SHZ1 The accuracy 
of the numerical values is high enough so as to render the 
corresponding curve indistinguishable from the 'correct' one 
and is plotted in its entirety, i.e. from the thin ring limit 
right up to the mass-shedding limit. The curves to first and 
third order were drawn by taking the expression for Qi, q , 
Cl 2 , M and J to first and third order respectively, inserting a 
numerical value for a and then taking the appropriate com- 
bination of these numbers. One finds in all three plots that 
the third order brings a marked improvement as compared to 
the first one, but that the behaviour near the mass-shedding 
limit is not particularly well represented. 
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M 

Qo 




Figure 6. The dimensionless mass divided by the outer radius 
is plotted versus the radius ratio for rings with polytropic index 
n = 1. This quantity tends to the value M/q q = 27rjoiJl(ioi) = 
7.84 ... in the thin ring limit, which is marked by a tick. 



20 



- - first order 

third order 

• • • • numerical 




Figure 7. The dimensionless angular momentum divided by the 
square of the outer radius is plotted versus the radius ratio for 
rings with polytropic index n = 1. This quantity as a function of 
radius ratio tends logarithmically to infinity. 



6 SOLUTION FOR AN ARBITRARY 
POLYTROPIC INDEX 

As was mentioned above, our generalized Lane-Emden equa- 
tion for /ioo can only be solved in closed-form for n = 1. For 
other polytropic indices, the iterative method presented here 
was applied with the help of numerics. By describing the 
unknown density terms fiik by Chebyshev polynomials and 
expanding all the quantities involved in terms of A, equa- 
tions can be formulated for purely numerical coefficients. 
The equations of the approximation scheme described in 
Section [5] must be fulfilled, whereby the ODEs for p,^ are 
evaluated at collocation points of the Chebyshev polynomi- 
als. In general, the density functions /ioo are not analytic at 
x — a, meaning that high order polynomials may be neces- 
sary to find a good approximation of the function desired. 
We none the less chose this method, since the equations in- 
volve integrals over the density for which one end-point of 
integration contains the unknown surface function r s , mak- 
ing their polynomial representation particularly useful. 

If one is only interested in determining a, /3n and Q.2, 
then it is not necessary to combine such numerical and alge- 
braic techniques and one can choose any numerical method 



Table 1 . The values of expansion coefficients for the surface func- 
tion r s and for the squared angular velocity are provided up to 
first order for different polytropic indices n. The value of a for 
n = can be found by solving J39D with n = and the conditions 
gj/^OoU=0 = an d /ioo(0) = 1 and then locating the first zero of 



A»oo- 


n 


a 











0.5642 





A + 3/4 




0.5 


0.7566 


-0.03537 


0.6371A + 0.5575 




1 


0.9594 


-0.07708 


0.4318(A + 1) 




2 


1.427 


-0.1731 


0.2169A + 0.2711 




5 


3.750 


-0.5118 


0.03614A + 0.07228 




10 


15.18 


-1.126 


(2.401A + 7.804) x 10" 


-3 


20 


207.6 


-2.375 


(1.362A + 7.829) x 10" 


-5 


30 


2.661 x 10 3 


-3.625 


(8.487A + 70.02) x 10" 


-8 


40 


3.337 x 10 4 


-4.875 


(5.468A + 58.78) x 10~ 


10 


50 


4.142 X 10 5 


-6.125 


(3.577A + 47.40) x 10~ 


12 



for solving the ODEs. One begins by solving equation (|39|l 
numerically for the desired polytropic index n, prescribing 
the 'initial conditions' /ioo(0) = 1 and ^fioo\ x _ Q = 0. For 
spherical polytropic fluids, a surface of vanishing pressure is 
known to exist only for n < 5, where the surface for n = 5 
extends out to infinity. The situation for polytropic rings is 
quite different - it seems that arbitrary polytropic indices 
are possible! Numerical solutions to (|39[) indicate that the 
density function /too indeed falls to zero for large n. The 
value of x at the first zero of the solution is a. One then 
proceeds to solve equation 



dx 2 



+ 



ld_ 

x dx 



1 



x- 



/in + 



4im _ n _ 
n + 1 



i - 



1 d/2 o 
dx 



(86) 



for /in with the condition /in(0) = and where ^Mii| _ 
has to be chosen so as to fulfil the centre of mass condition 



((9| to first order, which reads 



0=1 Hnx 2 dx = n / fiQ 1 fi,±ix 2 dx 



= a 



d/ti 



da; 



a/tii 



(«> + §/ 



(87) 



/ioo x dx. 



The constant /3n can then be found using equation {19 
which now reads 



An (a) / d/ipo 
a \ dx 



Pu = 



and Q.2 is taken from (|32[) . The behaviour of these coeffi- 
cients as they depend on the polytropic index n can be found 
in Table [T] The table suggests that a — > oo and Q.2 — > ex- 
ponentially in n for n — » oo, which is indeed known to hold 
l|Ostrikerlll964al lbl). The behaviour of the specific kinetic en- 
ergy of a particle in the ring, proportional to a 2 Q,2 to leading 
order, will be discussed in the next subsection together with 
the behaviour of 0u for large n. 

Before doing so, we provide a comparison of precise nu- 
merical values for various physical quantities with their first 
order equivalents in Table [2] One can see that the accuracy 
of the method does not depend strongly on the polytropic 
index and that relative errors are within a few percent for 
rings with a radius ratio of 0.9. 
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Table 2. Physical quantities to the first order in <r (bold faced type) are compared to the correct, numerically 
determined values (normal type) for given polytropic index n and given radius ratio Q\/ Qo = 0.9. 



n 


M 


n 2 




J 


P 


T 


W 


0.5 


103 


0.0216 




3.14 X 10 3 


59.2 


231 


-640 


0.5 


105 


0.0213 




3.21 x 10 3 


59.9 


235 


-650 


1 


143 


0.0151 




5.84 x 10 3 


89.2 


359 


-986 


1 


144 


0.0150 




5.95 x 10 3 


90.1 


364 


-999 


3 


356 


4.56 x 10 


-3 


3.53 x 10 4 


264 


1.19 x 10 3 


-3.17 x 10 3 


3 


359 


4.52 x 10" 


3 


3.58 x 10 4 


266 


1.20 x 10 3 


-3.21 x 10 3 


5 


714 


1.58 x 10 


-3 


1.44 x 10 3 


570 


2.86 X 10 3 


-7.44 x 10 3 


5 


720 


1.56 x 10- 


3 


1.46 x 10 3 


575 


2.89 x 10 3 


-7.51 x 10 3 



7 THE LIMIT OF INFINITE POLYTROPIC 
INDEX 



As n tends to infinity, the polytropic equation (|34[) shows us 
that pressure and density are proportional 



(89) 



a case sometimes referred to as 'isothermal' because such 
an equation holds for an ideal gas at constant temperature. 
Inserting this into equation Q35p at leading order and again 
using the dimensionless coordinate x yields 



Attx^qo + — ( x ^ hi moo j = <>. 



(90) 



The solution to this equation with our normalization 
[ioo(x = 0) = 1 reads 

4 



Moo 



7T [X + 



(91) 



The density and pressure fall to zero as x — > oo <4> r — > 
co. Integrating over the density to calculate the normalized 
mass, one finds to leading order 



M = 4-7T b J fioox dx — 4-Kb, 
o 



(92) 



which can also be read off from equation (123 directly, by 
making use of M = P, which is self-evident upon taking 
(|89p into account. In Fig. [SJ the behaviour of M /4kI> can be 
followed from the homogeneous case, n = 0, right up to the 
isothermal limit n — > oo. 

Making use of (|9ip . we find that 



lim g 

n — >oo 



0, 



(93) 



where g was defined in (|25[) . It thus follows from (|31[) that 

n\ 1 



lim (p 

71 — >CO \ 



n + 



(94) 



as already suggested by the results of Table [T] We can then 
see that the specific kinetic energy a 0,2 tends to infinity 
such that for fixed A 



lim 

rwoo n + 4A 



= 1. 



(95) 



From the fact that |/3n| tends to infinity, we can con- 
clude that the range of a values for which the first order pro- 
vides a good approximation shrinks to the point a = 0. This 
provides us with evidence suggesting that the deviation in a 
ring's cross-section from a circle becomes more pronounced 



1.0 



0.9 - 



M_ 
4irb 



0.8 - 



0.7 



0.6 - 



0.5 



0.0 0.2 0.4 0.6 



0.8 



1.0 



n+1 



Figure 8. The dimensionless mass divided by 4irb in the thin ring 
limit is plotted versus n/(n+l) over the whole range of polytropic 
indices n £ [0, oo). The points for n = 0, n = 1 and n — > oo are 
known analytically and the remainder of the curve was generated 
by solving the equation for /ioo numerically and making use of 
equation (1206 . 



at a given radius ratio as n is increased. The value for Qi/ Q 
at which one reaches the mass-shedding limit presumably 
tends to 1 as n tends to infinity. 
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IPrudnikov et al.lll990l ) 

oFi (&,-*) = V(b)z (1 - b)/2 J b ^{2^) 

The above double integral yields 
i i 

J (2^/ht^) zJi{2y/t^z) 



2y/tlt 2 (l -ti) 2^Jz(l-t 1 ) 

1 f Jq (2^771) ^ 2 
t 1 (l~t 1 ) dtl = J °^ ) 



dti 



thereby proving (|A1|I . 



APPENDIX A: AN IDENTITY RELATING 
HYPERGEOMETRIC TO BESSEL FUNCTIONS 

In order to prove the identity (|84p . we prove the more general 
identity 



40 



5 5 



2' 2' 



77^2^3 77, -;3,3,-;-z -7:2^3 ~, -;2,2, 



3 3 



2' 2' 



(Al) 



for an arbitrary complex number z, from which (|84|l follows 
immediately. 

We begin by using the differentiation properties of the 
hyper geometric functions, e.g. 7.2.3.47 in IPrudnikov et al.l 
(1990), to write 



40 



2F3 



5 5 



7 



, -; 3, 3, -; — z ) z — - 2F3 



2' 2 



_d_ 

dz 



2F3 



2 

1 I-l 1 l.-z 
2'2' ' '2' 



2'2' ' '2' 



2F3 



5 3 -22 l--z 

2'2' ' '2' 



(A2) 



With the integral identity 7.2.3.11, the term to be differen- 
tiated can be written as 



2F3 



i i-l 1 l-z 

2'2' ' '2' 



2F3 



3 - *-2 2 5 -- Z 
2'2' ' '2' 



1 ff Jo (2^777^) _ z Ji (2^/77^) 



(A3) 



dtj dt 2 



where we have made use of the identity (e.g. 7.13.1.1 in 



